Spin noise of itinerant fermions 
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. We develop a theory of spin noise spectroscopy of itinerant, noninteracting, spin-carrying fermions 

in different regimes of temperature and disorder. We use kinetic equations for the density matrix 
in spin variables. We find a general result with a clear physical interpretation, and discuss its 
dependence on temperature, the size of the system, and applied magnetic field. We consider two 
classes of experimental probes: 1. electron-spin-resonance (ESR)-type measurements, in which the 
probe response to a uniform magnetization increases linearly with the volume sampled, and 2. optical 
Kerr/Faraday rotation- type measurements, in which the probe response to a uniform magnetization 
increases linearly with the length of the light propagation in the sample, but is independent of the 
' cross section of the light beam. Our theory provides a framework for interpreting recent experiments 

, on atomic gases and conduction electrons in semiconductors and provides a baseline for identifying 

the effects of interactions on spin noise spectroscopy. 
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I. INTRODUCTION 



Currently, there is much interest in studying the physics of nano-scale structures. In measurements of a response 
function by pump-probe experiments on systems of decreasing size, the signal decreases more rapidly than the noise, 
and thus the signal-to-noise ratio decreases with decreasing system size. The fluctuation-dissipation theorem, which 
relates a response function to a correlation function obtained from noise measurements, enables us to change this 
d • problem into a useful tool. An additional advantage of noise measurements is that they often disturb the system less 

than experiments that measure the response of the system to an external perturbation. 
I , There have been a number of experiments studying spin properties of systems using spin noise. Spin noise has been 
measured in systems of spins whose position is fixed in space: atomic nuclei, spin glasses^, magnetization modes 
in magnetoresistive heads^, and electrons and holes in self- assembled quantum dots 4 . There have also been recent 
measurements of spin noise of itinerant spins in hot atomic gases^i, cold atomic gases ? ' in rt-doped bulk GaA o 10 ' 11 , 
and in n-doped GaAs quantum wells 1 ^. Localized spin noise measurements on nanostructured systems using STM 
techniques have been discusse d 13 ' 14 . The experimental setup of spin-noise spectroscopy in semiconductors has been 
optimized in Refj 1 ^ and has been used to measure spatially resolved doping concentration in GaAs 1 ^. Motivated 
by these experiments with itinerant spins, we have developed a theory of spin noise of itinerant fermions in different 
regimes of temperature (degenerate/classical statistics) and disorder (ballistic/diffusive motion). We consider the case 
, of noninteracting particles as a benchmark for comparison to experiments, so that we can then identify the effects 
■^j- ' of interactions on spin noise. We find a general result that holds in the different regimes, which has a clear physical 
interpretation, and we show how it follows from kinetic equations for density matrix in spin variables. We consider 
I ■ two classes of experimental probes: 1. electron-spin- resonance (ESR)-type measurements, in which the probe response 
° to a uniform magnetization increases linearly with the volume sampled, and 2. optical Kerr/Faraday rotation-type 
. measurements, in which the probe response to a uniform magnetization increases linearly with the length of the light 
^ ' propagation in the sample, but is independent of the cross section of the light beam. 
• i— i . The outline of the paper is as follows: In Sec. HH we present the general result and its interpretation. In Sec. IIIIl we 
show how the noise spectrum behaves as a function of temperature, system size, and magnetic field. In Sec. IIV1 we 
present the derivation of the result from kinetic equations with details of calculations in some limiting cases presented 
in the Appendix B. In Sec. |Vl we present our Conclusions. The spin noise power spectrum of spin J, which motivates 
our general result, is derived in the Appendix A. 



II. GENERAL RESULT 



In Fig. [T] we show the generic setup of the experiments considered. We study the noise of spin magnetization in 
the z direction in the presence of an applied constant magnetic field in the x direction that splits the two spin energy 
levels by the Larmor frequency u)l. The purpose of the constant magnetic field is to shift the noise spectrum away 
from zero frequency. Noise measurements close to zero frequency are difficult because of the presence of ubiquitous 
1// noise. The magnetic field is chosen so that ldl is larger than the linewidth of the spin noise spectrum. For charged 
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FIG. 1: Schematic view of the experimental setup. The experiment measures the noise of spin magnetization in the z direction 
in the presence of an applied constant magnetic field in the x direction that splits the two spin energy levels by the Larmor 
frequency lul- The system has thickness L in the z direction, and is extended in the x — y direction. Noise in the part of the 
system with transverse size 7? and cross section A ~ R 2 is probed. 



fermions, we neglect coupling of magnetic field to the orbital motion, that is, we consider the case when the cyclotron 
orbit is longer than the dimensions of the probed region. The system has thickness L in the z direction, which in the 
optical experiments is the direction of light propagation. The system is extended in the x — y direction. Noise in the 
part of the system with transverse size R and cross section A <~ R 2 is probed. For a uniformly spin-polarized sample, 
the magnitude of the probe response in bulk measurements such as ESR scales as the volume of the probed region, 
that is, Signal oc R 2 L. By contrast, for a uniformly spin-polarized sample, the magnitude of the probe response in 
optical Kerr/Faraday rotation measurements scales as the thickness L, but is independent of the cross-sectional area 
A, that is, Signal oc L. 

In general, the spin noise experiment will measure the noise power spectrum of a quantity Q proportional to the 
instantaneous electron spin polarization 

Q(t)=CM z (t)(GF), (1) 

where M z (t) is the operator of the z— component of the instantaneous electron spin polarization in the probed volume 
at time t, related to the spin-density operator s z (r,t) by 



M z (t)= J d A rs z (v,t), (2) 

AxL 

C is a fixed coupling constant, and (GF) is a geometric factor. For bulk measurement such as ESR, the geometric 
factor is unity. For optical Kerr/Faraday rotation measurements, the geometric factor is 1/A. Hence, in either case, 
in order to calculate the spin noise power spectrum, we need the Fourier transform of the correlation function 

S zz (t 2 -h) = ^({M z (t 2 ),M z (h)}), (3) 

where Here, {, } denotes the anticommutator, and () is the equilibrium ensemble average. 

In equilibrium, the fluctuation-dissipation theorem relates the noise power to the imaginary part of the corresponding 
susceptibility 

S zz M = -coth^ J AAxL^-r^)- -coth^Ax L^ z (q = 1/ R;uj), (4) 

AxL 

where 



X zz (r 2 - ri;t 2 - *i) = -id(t 2 - t 1 ){[s z (r 2 ,t 2 ), s z (r 1 ,t 1 )}). 



(5) 
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FIG. 2: General qualitative form of the noise power spectrum. It is peaked at the Larmor frequency lul- Its width is 
approximately equal to the larger of the inverse travel time t^ 1 and the inverse spin- flip time r s _1 . The height is given by the 
magnetization Mx multiplied by the thermal factor cothLUL/2T and divided by the width. 



Here, \zz is the spin susceptibility. 

We consider different regimes of temperature and particle motion. The calculations for the various cases are given 
in Sec. IIVI The results of the calculation in all these cases have the following general form 



S zz (uj) ~ ^-coth 1 



M 



(0) 



2""" 2Tg(t^,r s - l y 



(6) 



Here, is the equilibrium magnetization of the probed region caused by the constant magnetic field in the x 

direction, t tr is the travel time it takes the spin-carrying fermion to move across the probed region (distance i?), t s 
is the spin-flip time, / is a function of unit height and width, peaked at zero, and g is a function whose value is 
approximately equal to the greater of the two arguments. The detailed form of the functions / and g depends on 
temperature and disorder. The form of the functions / and g for the various cases will be discussed in section ITVl 
Equation [5] is the main result of this paper. Its form agrees with the spin noise power spectrum of a single spin J-see 
Appendix A. 

In Eq. [5J for ballistic transport, 



R/v, 



(7) 



where v is the Fermi velocity y / 2EF/m in the degenerate regime T « Ef or the thermal velocity y / 2T/m for 
T » Ep. For diffusive transport, 



ttr = R 2 /D, 



(8) 



where D is the diffusion constant. 

In Fig. [21 we show the schematic behavior of the spin noise power spectrum S zz {oj). It is peaked at the Larmor 
frequency lol- Its width is approximately equal to the larger of the inverse travel time t^ r and the inverse spin- flip 
time t.T 1 . The height is given by the magnetization M% multiplied by the thermal factor cothwi/2T and divided 
by the width. 

The result j6| has a simple physical interpretation. The scale of the response in the z direction is set by the initial 
polarization in the x direction. The width is given by the inverse of the time to lose spin coherence, either due to a 
spin flip or by moving out of the probed region. This inverse time divides the magnetization to give the noise power 
the dimension of inverse frequency. The peak in the noise power spectrum is centered at the Larmor frequency tJi,. 
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III. IMPLICATIONS 

We discuss consequences of Eq. (|6]). Because Mx grows linearly with the probed volume, and the other quantities 
in Eq. © are independent of L, the height of the power spectrum grows linearly with L while its width is independent 
of L. 

The dependence of the noise power spectrum on R is different for the ballistic and diffusive motion. In the ballistic 
case, the height of the noise power spectrum behaves as 

R 2 

S zz (uj l )(x— — — , (9) 



and its width behaves as 



gitt^T-^K^+T- 1 . (10) 

Thus, for R < vt s , the height of the noise power spectrum scales as i? 3 , and its width as 1/R, whereas for R > vt s , 
the height of the noise power spectrum scales as R 2 and its width is independent of R. In the diffusive regime, the 
height of the noise power spectrum behaves as 

R 2 

S zz (uj l ) oc D ■ __ 1 , (11) 



W i T s 



and its width behaves as 



9(t7r\rr 1 )<x^+T- 1 . (12) 

Thus, for R < \/Dt Si the height of the noise power spectrum scales as R 4 and its width as 1/R 2 , whereas for 
R > \/Dt s , the height of the noise power spectrum scales as R 2 and its width is independent of R. For both the 
ballistic and diffusive transport in the long-i? limit, we recover the behavior of static spins flipping on the time scale 
t s . For optical Kerr/Faraday rotation experiments, there is an additional factor 1/A 2 in the height of the noise power 
spectrum for the geometrical factor, Eq. [1] The scaling behavior is summarized in Fig. [3] The scaling crossover is 
most pronounced in optical measurements of the noise in the diffusive regime. 

We consider the dependence of the noise power spectrum on the external magnetic field, that is, on the Larmor 
frequency u>l- The width is independent of ui^. Hence, instead of studying the width and the height of the noise 
power spectrum separately, we study the integrated spin noise power spectrum 

oo 

S zz ee / dujS zz (uj), (13) 



n 



which scales as the product of the width and the height. ^From Eq. [SJ we see that 

^excoth^MW, (14) 

that is, it is independent of the characteristic rate g{t^ , t" 1 ). In the classical regime, the magnetization M x ^ depends 
on ljl as 

MWcctanh^. (15) 

In the degenerate regime, m| ^ also has a linear dependence on u>i for small lol, cut off at the Fermi energy Ep. We, 
therefore, approximate the field dependence of by 

M<°> cxtanh ^r^, (16) 

max(2i , hjp) 

giving 

S 7 , oc coth — - tanh , L „ ; . (17) 

2T m-Ax(2T,E F ) V ' 

In the classical regime, 2T/Ep > 1, the ^-dependence of the two factors in ([17]) cancels, and S zz is independent 
of the magnetic field, in agreement with the measurements of spin noise in Rb vapors!. In the degenerate regime, 
2T/Ep < X, there is an intermediate regime of magnetic fields 2T < lol < Ep, where S zz grows linearly with field, see 
Fig. HI The spin noise measurements in GaA a 10 i n ' 12 were done at temperature T >> lol, so S zz is field independent. 
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FIG. 3: Scaling of the height (solid line) and width (dashed line) of the noise power spectrum as a function of R in the ESR 
measurements (first row) and optical measurements (second row). The first row are the log- log plots of formulas @, (|1U[1 and 
(If f I) . (|f 2|) ; the second row are the log- log plots of the same formulas with the height divided by R 4 . The heights and widths 
are normalized to their value at R — vt 3 for the ballistic case (the first column) and at R — y/ Dt s for the diffusive case (the 
second column). The scaling in the different regions is denoted at each curve. The crossover is most pronounced for the optical 
measurement in the diffusive regime. 
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FIG. 4: Solid line: the magnetic-field dependence of the integrated noise power spectrum S zz m the classical (2T/Ef = 3, 
left panel) and degenerate (2T/Ef = 0.3, right panel) regime. Dashed lines: the factors coth ^ (upper dashed lines) and 
^ ann max(2T Bp) 1 ^ ne vame °f ^ s normalized to its value in the high-field limit. In the classical regime, the field dependence 
of the two factors cancels, and S zz is field independent. In the degenerate regime, S zz grows linearly with field in the region 
2T < lul < E F . 
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IV. CALCULATION 



We now turn to a detailed justification of the above results. We calculate the susceptibility \zz as a linear response 
of the spin density (s z (r, t)) to an external potential <p(r, t). The Hamiltonian describing the coupling is 



Ha, 



d 3 r0(r,t)s z {r,t). 



(18) 



In terms of the electron field operator in the Heisenberg representation, ip(r,t), the spin density operator s z (r,t) is 
given as 



s z {v,t) = ^(r,t)^(r,f), 



(19) 



where a z is a Pauli matrix. In order to calculate the linear response, we construct and solve the kinetic equation for 
the density matrix in the Wigner representation 



AT«P' f Unr--,fW; a r+-,t 



p a/3 (p,r, t) 

In terms of this density matrix, the spin density is 

(«»(*,*)) = 

The Hamiltonian consists of three terms: a non-interacting term 

#o = 



^tr(p(p,r,^). 



d 3 r (^-W 0, W(r, t) - u l ^(t, t)^(v, t) ) 
\ 2m 2 / 



(20) 



(21) 



(22) 



a coupling term to the external potential (|18[) . and a scattering term, which determines whether the particle motion 
is ballistic or diffusive. The first term on the right-hand side of Eq. [22] is the kinetic energy, and the second term 
describes the interaction with the applied magnetic field. The kinetic equation for the density matrix has the following 
form 



d tPa p{p,v,t) - i / d 3 fe-* p f 



H a + H^^\r~ -,t)ip a \ r+-,t 



where I s describes the effect of scattering. Substituting (fT5)) and (f2"2"j) into (f2"3")l . we find 



= 1, 



(23) 



IUJ L 



( dt+ £ . V ) p(p,r,t) 
To obtain linear response, we write 



y,P(P,r,t) 



r - At p(p,r,i)— - cj> r+-d p ,t —p(p,r,t) = I 



p(p,v,t)=p^{p)+6p(p,r,t). 
Here, p^°'(j>) is the equilibrium density matrix corresponding to the Hamiltonian Ho, that is, 

^ 0) (P) = \ (n (ik9) - f ) + n (m + l f))-o + \(n (f(p) - f ) - n (f(p) + f )) <r. 



(24) 
(25) 



(0)/ 



Here 



and 



"(0 = 



e r +1 



e(p) 



2m 



(26) 



(27) 



(28) 
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The particular deviations from equlibrium we consider have the form 

8p(p, r, t) = 8p y (p, r, t)o y + 5p z (p, r, t)a z = <5p(p, r, t) ■ er. (29) 

The scattering contribution to the kinetic equation is different for the ballistic and diffusive motion. We first 
consider the case of the ballistic motion, in which 

I s = --Sp(p,r,t). (30) 

ft describes spin relaxation with relaxation time t s . 

To obtain the susceptibility at positive frequencies, we write the kinetic equation for the circular component of the 
density matrix 

Sp+(jp, q, w) = J d 3 rdte- i<l - r+iut {8p z - i8p y )(p, r, t). (31) 

The equation for the opposite circular component 5p- gives the susceptibility at negative frequencies, which gives an 
equivalent result. The kinetic equation becomes 

-i w- q+— 5p + (p,q,w) = i<t>(q,w) ^—z 1 — (32) 



m r, _ 
To obtain the susceptibility, we solve for 

8p+(q,u) = J ^<5/)+(p 1 q,w) ) (33) 

which gives 

A{ „ u > o) -\ f ^ " « <P - » - V) - » ft + I) + *) . (34) 

4 7 (27T) d w _ WL _IL.q + JL. 

The superscript i? denotes the case of ballistic motion. 

For diffusive motion, the fermion momentum relaxes on a rapid time scale r << r s . In this case, the scattering 
contribution to the kinetic equation is 

Is = -—Sp(p,r,t) - -(p(p,r,t) - p eq (p;p(r,t))), (35) 

where p e q{p] p(r, t)) is the momentum equilibrium distribution consistent with the local number-spin density p(r, t), 

d 3 p 



(27T)' 



iPeq{p;p{r,t)) = p{r,t). (36) 



For <5p(p, r, i) given by Eq. 



1 , , . I, f .1 n ($(p) - - n (f(p) 



To obtain the susceptibility at positive frequencies, we write the kinetic equation for the circular component of the 
density matrix 

.( P _< t\_ . . ... , P^ 0) (P + | ) - pf_ (P - | ) - pi 0) (P + | ) - pl 0) (P - |) 

-t q+ — + -Up+(p,q,w) = i0(q,w) * ^ * ^— * ^ * 

\ m t s t J 2 

+ ; /^(»(« P )-^)-..(«P) + *)) P+(q ' (38) 

To obtain the susceptibility, we solve for 

d 3 p 



8p + (q,w) = J ^-^Sp + (p,q,u) 



(39) 
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R < V S T 

(ballistic) 


R > V S T 

(diffusive) 


T < E F 
(degenerate) 


Trapezoidal (for ojl/2 < Ef ) 
Parabolic (for wl/2 > Ef) 


Lorentzian 


T > Ef 
(classical) 


Gaussian 


Lorentzian 



TABLE I: Summary of the shapes of the noise spectral lines in the four regimes in the limit t 3 1 — > 0. 



which to lowest order in r gives 

>o) 



4 



(£(p) 



2 . 



»(e(p) + ^)) 



<jJ - OJL + i{r s 1 + D{uJL)q 2 ) 
The superscript D denotes the case of diffusive motion. Here, 

/AiMe(p)-f)-"(c(p)+f)) 



(40) 



(41) 



In the limit ojl — > 0, D{uj£) is equal to the usual diffusion constant. For the classical distribution n(£) in the high- 
temperature limit, D(ujl) is independent of u>l for any value of ojl. At T = 0, the high-field limit D(u;l/2 > Ef) is 
reduced compared to the low-field limit D(u>l — > 0) by a factor of 2 2 / 3 x 3/5 = 0.95. Thus we consider D independent 
of the magnetic field, and drop the col dependence in the following. 
Substituting Eq. [34] into Eq. |4] gives 



^M~|coth^AxL 



d 3 p 
(2tt) 3 



(Hp 



2R) 



¥)-"(C(p- 



2R> 



2 , 



1 

7TT S 



(42) 



- — • 



for the case of the ballistic motion. Here, e is the unit vector in an arbitrary direction. Substituting Eq. 0D]into Eq. 
S] gives 



S£(w) ~ - coth —A x L 



d 3 p 



(e(p)-^)-n(e( P ) + ^) 



i^-i+D/i? 2 ) 



( W -^ L )2 + (r7 i +^/i? 2 ) 2 



(43) 



for the case of the diffusive motion. Evaluating the integrals in Equations [42l and l43l see Appendix B, gives results 
of the form of Eq. [6J with specific forms for function f Uw — WL)/g(f^ r , t" 1 )). The forms of the function / in the 
various cases are summarized in Table I. 

In the diffusive case, the function / is a Lorentzian for all temperatures, and 



g(t tr \T a x ) = C + 



(44) 



The function / in the ballistic case (|34[) is more complicated, so we consider various limiting cases. In the limit 
T.T 1 >> v/R, where v is the Fermi velocity at T = and the thermal velocity \J2T jm at T > Ep, we can drop 
p/m • q from the denominator of (|34p . so the numerator becomes q independent, and we obtain a Lorentzian with 
the width T.T 1 , like in the diffusive regime for t^ 1 » D/R 2 . In the limit of r^ 1 dominating the inverse travel 
time through the probed region, the spin noise cannot distinguish between the ballistic motion and diffusive motion. 
We can distinguish between these cases in the opposite limit r s _1 — * 0. Then the line in the diffusive regime is still 
a Lorentzian, but now with the width D/R 2 . In the ballistic case, at T = 0, the line has a trapezoidal shape for 
ujl << Ef and a parabolic shape for ljl >> Ef- In the classical limit T >> Ef, the line has a Gaussian shape. 



V. CONCLUSIONS 



Motivated by spin noise spectroscopy measurements, we developed a theory of spin noise of itinerant fermions 
in different regimes of temperature and disorder. We found a general result with a clear physical interpretation, 
and showed how it follows from spin kinetic equations. Our theory provides a framework for interpreting recent 
experiments on atomic gases and conduction electrons in semiconductors and provides a baseline for identifying the 
effects of interactions on spin noise spectroscopy. 
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Appendix A 

In this Appendix, we derive the spin noise power spectrum S zz (u>) for a single spin J in the magnetic field applied 
in the x direction: 

s zz (t 2 -h) = \({J z (h),Jz(h)}) 

= i( e -iu, £ (t 2 -t 1 ) + e i Wi (i 2 - 4l )^| J+) j_yj (45) 

8 

where 

J± = J z TUy (46) 
The thermal expectation value of the anticommutator is related to the thermal expectation value of the commutator 

({J + ,J_})^coth^([J + ,J_]>. (47) 

Using the commutation relation 

[J+,J-)=2J X , (48) 

we arrive at the result 

S„(w>0) = |coth^(J B )%-w L ). (49) 



We see that Eq. 032 has the form of Eq. [5] with 

1 

steW 1 )" V</( i;, L -r: '■) 

The formula for magnetization 



^—rj( u-i^-is )=6(u>-u> L ). (50! 



(J x ) = ¥±1 coth (2J + 1 ^ - i coth ^ (51) 

implies that in the special case of J = 1/2, the temperature and field dependence cancels between cothwi/2T and 
(J x ) giving 

S„(u>0) = j8{u-u L ). (52) 
Appendix B 

We present the details of calculations for the ballistic case with r" 1 — ► in three different limiting cases. 

A. T = 0, ul « E F 
In this case, we can linearize around the Fermi surface. From f|34|) . we now get 



l 



v F qc + uj L \ ( v F qc + uj L 

-t, ^ ) ) d{u;- {VFqc + ujL)) 



7r N uj l uj 6(v F q U) L \) ,_„, 

2 X 2 X uj l 2v F q ' 
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We thus obtain ([6]) with the magnetization proportional to the magnetic field via the Pauli susceptibility, 
N lol/2, ballistic transport time t tr — R/vp, and / with a trapezoidal shape. 



B. T = 0, lu l /2 > E F 



In this case, the fermions are fully polarized, so 

oo 1 



u; — ujl qc 

m 



1 / p 

- W L 

2 V m 



p \ p 2 + {q/2f 
qc) - 



2m 



! 

1 m 2 
167T q 







where we neglected q 2 /2m compared to the Fermi energy /x + u>l/2. Using the formula for particle density 

1 „ 



with 



we find 



7T 71 3 1 



Xb(<1,u > 0) = -- x - x - 



2 2 Aqp F /m 



1 - 



qp F /m 



(54) 



(55) 



(56) 



(57) 



The magnetization is now saturated at one half times the number of the fermions in the probed region. The sus- 
ceptibility is non-zero because it is transverse to the external magnetic field [zz response function with the external 
magnetic field applied in the x direction). The transport time is still equal to R/vp. The line shape / is now parabolic. 



C. T> E F 



In this case, we approximate 



with \i < 0, so © gives 



Xb(9,^>0) 



n(0 = e-^-^'T 



(58) 



J p 2 dp J dee ' sinh ( vqc ^ T " L \ 5{uj — (pqc/m — uj l )) 

ir o -l v ' 



J p 2 dpe cosh ^ 
o 



tt n u L e 2T ^/<- 
- — x — tanh — x 
2 2 



(59) 



2T ^2nTq 2 /m 

Thus, the equilibrium magnetization is n/2 tanhcji/2r, the transport time is R/v, where v is the thermal velocity 
y/2T/m, and the line shape / is Gaussian. 



* Electronic address: simonkos@kfy.zcu.cz 
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